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Abstract 

We have calculated the temperature dependence of the frictional drag be- 
tween spatially separated quantum wells with parallel two-dimensional elec- 
tron gases due to interlayer electron-electron interaction mediated by virtual 
exchange of acoustic phonons due to piezoelectric and deformation potential 
interaction. It is shown that the frictional drag is dominated by the piezo- 
electric coupling. According to our calculations the temperature dependence 
of the drag scattering rate divided by T 2 exhibits a pronounced peak which 
for the experimental situation and in agreement with the finding of T. J. 
Gramila, et al. Phys. Rev. B 47, (1993) 12957 is obtained at about T r* 2.1 
K. We ascribe the appearance of this peak to a change from small to large 

angle scattering in the virtual phonon exchange. 
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I. INTRODUCTION 



Double layer two-dimensional (2D) carrier systems as realized in semiconductor nanos- 
tructures are of interest with respect to fundamental phenomena as fractional quantum Hall 
effect and superconductivity because of the interlayer carrier interactioni. Depending on 
the width of the potential barrier that separates the two layers one has to distinguish the 
regime of quantum-mechanical coupling via tunneling for small barrier widths from that of 
interlayer coupling mediated by (direct or effective) electron-electron interaction for suffi- 
ciently large barriers. In the former case, the frictional aspects of super conduct ivityi'i and 
of the fractional quantum Hall effect&i have been studied. In the latter case, when the 
large barrier allows independent control of electron conduction in the two layers, the fric- 
tional drag effect manifests itself when a current driven along one layer (layer 1) induces, 
via momentum transfer, a drag voltage in the second layer (layer 2) under conditions that 
no current flows in this welliQ. The drag effect has been studied experimentally for two 2D 
electron systems (2DES), two 2D hole systems (2DHS), but also in situations with layer 1 2 
being a 2DES or 2DHS, respectively^^. Recently the effect of a magnetic field perpendic- 
ular to the layers has been investigated!^!. Theoretical work on the frictional drag effect 
was devoted to elucidating the microscopic mechanism of the momentum transfeJil~0. The 
calculated temperature dependence of the frictional drag due to direct Coulomb electron- 
electron interaction shows strong dependence on the inter-layer separation AlHii while in 
the experiment by Gramila et al.El, the temperature dependence of the observed frictional 
drag for separations at least up to A ~ 50 nm exhibits almost no dependence on A, and 
has been interpreted as being due to the exchange of virtual phonons. The frictional drag 



accomplished by the exchange of virtual phonons has been studied theoretically in |24j. This 
work, however, because of its conceptual complexity requires involved numerical computa- 
tions and does not lead to analytical expressions for the dependence of the frictional drag 
on the system parameters. Therefore, we present here an analytical derivation of the drag 
scattering rate for the particular case of asymmetric two- layer systems {n\ ^ n 2 ) coupled 
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by virtual exchange of acoustic phonons including its temperature dependence. We thus 
model the experimental situation investigated by T. J. Gramila, et al.El and find by using 
the corresponding system parameters the characteristic T-dependence of the drag scattering 
rate in quantitative agreement with these experimental data. 



II. THEORETICAL CONCEPT 

Our model as depicted in Fig. [I] consists of two parallel electron sheets separated by a 
distance A with electron densities n\ and n 2 and the layer extensions d\ and d 2 . According 
to experimental observations, a current l\ driven along the layer 1 by applying a voltage V\ 
leads by frictional drag to a voltage V 2 in the layer 2. The drag scattering rate in a double 
layered system can be defined as0 

— = 4- a) 

1~Drag 

where Ei, E 2 are the electric fields connected with the voltages V\, V 2 , and T\ is the transport 
relaxation time which determines via the mobility fi = eri/m* the current in the layer 1; m* 
is the electron effective mass. Following p3l we obtain by using the linearized Boltzmann 



equation for the drag scattering rate 
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where L is the normalization length, k the electron in-plane momentum, a the electron spin, 
q± = k — k' the transferred momentum, / the Fermi distribution function determined by 
temperature T and Fermi energy €p = fi 2 k 2 F /2m* with Fermi wave vector related to the 
areal electron concentration n, by k 2 F = 27m. W^ 2 ^v 2' ^ s the transition probability of two 
electrons from the states |1 >= \a±, k\ > and |2 >= |cr 2 , k 2 > with energies E\ and e 2 , 
respectively, into the states \V >= |ov> k\> > and \2' >= \a 2 ', k 2 > > with energies Ey and 
e 2 i. The index T refers to the type of effective electron-electron interaction which causes 
this transition. 



In this work we consider effective electron-electron interaction mediated by the exchange 
of virtual acoustic phonons related to the piezoelectric (T = PA) and deformation (T = 
DA) potential coupling. Using Fermi's Golden Rule, W^-^v 2 1 can be written in the form 
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where T 1 T 2 _ +1 , 2 , is the transition matrix element. Phonon mediated electron-electron in- 
teraction appears in second order perturbation theory with respect to the electron-phonon 
coupling. The transition matrix element of this process is visualized by the diagram in 
Fig. 0. For intrasubband scattering the vertex part of this diagram corresponds to 

T r (q) = Jm(q) / dzp(z) exp(iq z z) (4) 

with the subband electron density p(z) and q = (q±, q z ). B r (q) depends only on the absolute 
value of q 

B p \q) = B?\-\ = ^f, (5) 
B DA {q) = B?\ Bg A = g. (6) 

Here we describe electron-phonon interaction in the framework of the isotropic model with 
s being the longitudinal sound velocity, g the crystal mass density, eft and S are the piezo- 
electric and deformation potential constants, respectively, averaged over the directions of 
the vector q and the phonon polarizationsil'il. We assume also that all elastic parameters 
of the sample are the same so that phonons are not reflected from the interfaces separating 
different materials (for instance between GaAs and AlGaAs). The phonon propagator is 

with hu = Ei — 6\r, u q = sq and a summation over phonon momenta q corresponds to 
the internal phonon line (the dashed line in the diagram). In the above expression of the 
phonon propagator, r q is the lifetime of the intermediate phonon states with respect to all 
processes which are destroying these intermediate states. For virtual intermediate states 



there is no energy conservation at the vertices of the diagram, i.e. E\ + 62 = £y + £2' can 
be different from tkv and we may assume infinite lifetime r q for the intermediate phonons 
when calculating the drag rate. Then all integrals with respect to the energy parameters 
in Eq. (|2[) are to be understood in the sense of the principal value. Taking the summation 
over intermediate momenta q, we obtain for the squared-modulus of the transition matrix 
elements for the two considered mechanisms 
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where 
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and a is defined generally as a(q±,cu + i£) = (g 2 — (a; + i£) 2 /s 2 ) 1//2 where £ is the positive 
infinitesimal, and the branch cut for the square root is assumed to lie along the negative 
real axis. It is easy to check that with such definition of a, the expressions of the form 
factors ( |10|) and (|ll|) are valid both for q\ — uo 2 / s 2 > and q\ — uo 2 / s 2 < 0. In calculating 
the form factors I PA and I DA we have assumed that electrons are localized in symmetric 
infinitely high quantum wells with the width d, thus having pi(z) = (2/d) sm(irz/d) 2 and 
P2{z) = pi(z + A) as explicit forms for the electron density functions in layers 1 and 2, 
respectively. In Eqs. @ and (H) we introduce also nominal scattering times 
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with the numerical values obtained for a GaAs quantum well with the electron concentration 
n x = n 2 = n » 1.5 • 10 11 cm" 2 as in 0. Substituting Eqs. (|),(§), and (@) into Ec i- (0) and 
taking integrals over ky and over directions of the vector k 2 by exploiting ^-functions, we 
can represent the drag scattering rate in the form 
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In this equation e q± = h 2 q 2 _/(2m*) and the limiting energy e (q±) = (hu — e g± ) 2 / (Ae g± ) is 
obtained from the energy and momentum conservation laws. Now making use of the identity 

tt wi t( 1 ^ f( x ) - f( x + y) 

f(x)(l-f(x + y))= (16) 

and taking integrals in Eq. fll5|) with respect to the energy parameters, we reduce the 
expression for J-" T (T) to a one-dimensional integral of the form 
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where PolyLog[i>, z] gives the polylogarithm function of the order 
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PolyLog[u, z] = Li u (z) = — - / rft (18) 
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and 
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For the special GaAs system with n ~ 1.5 • 10 11 cm -2 we have 7 ~ (7.6 K) jT and :rp ~ (62.4 

K)/T. One can see from Eq. ( |TTD that in the framework of the adopted approximation when 

the lifetime of the intermediate phonon states is infinite, i.e. r^ 1 = 0, the drag scattering 

rate, Eq. (|14|), is independent of the inter-layer spacing A. This means that our theory refers 
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to the experimental situation when A is smaller than the phonon mean free path, Ao = sr q , 
associated with the finite phonon lifetime. For A > A , we have to consider the finite r q in 
the phonon propagator, which would lead to an exponential dependence of the scattering 
probability of W oc exp(— A/A ). 

III. RESULTS AND DISCUSSION 

Using Eqs. ( |14D and fllTD we plot the temperature dependence of the drag scattering rate 
divided by T 2 for identical layers in a symmetric GaAs/AlGaAs double quantum well struc- 
ture as in [l^ coupled both by PA- and DA-phonon mediated electron-electron interaction 
(Fig. |3|) . It is seen from this figure that the contribution of the DA interaction is smaller 
(approximately by one order of magnitude) than that of the PA interaction. This is because 
the scattering probability due to long-range piezoelectric interaction includes an additional 
factor q~ A with respect to the short-range deformation potential interaction. This difference 
is strongly pronounced at low temperatures^ where the Pauli exclusion principle restricts 
electron-phonon scattering processes in a Fermi gas to those with small angles. 

Our calculations correspond to the experimental situation of Ref. |l^ where the drag rate 
was measured for samples with the inter-layer spacing A = 17.5, 22.5 , and 50 nm with the 
electron concentration n — 1.5 • 10 11 cm -2 in a temperature range from 1 to 7 K. Starting 
from the lower edge of this interval, where the drag rate is small because it is dominated by 
small angle scattering connected with a small momentum transfer, the drag rate increases 
with increasing temperature. Towards the upper edge of the interval, large angle scattering 
becomes dominant which reduces the increase of the drag scattering rate in temperature. 
This gives rise to a peak of the function 1/ (r^ rag T 2 ) both for drag caused by the exchange 
of DA- and PA-phonons at approximately the same temperatures as shown in Fig. |3[ The 
maximum drag rate of about 1 s _1 K~ 2 due to combined effect of DA- and PA-phonons is 
obtained at a temperature near T p « 2.1 K in good agreement with the experimental results 
for the samples with the inter-layer spacing A = 17.5, 22.5, 50 nm0. Our theory is not 
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applicable for much larger interlayer spacing (as e.g. for the sample with A = 500 nm in Ref. 
[13]). We believe that in this sample due to other scattering mechanisms, the intermediate 
virtual phonon states which effectively realize electron-electron interaction, are destroyed 
and this leads to the exponential suppression of the drag scattering rate observed in the 
experiment. 

In Fig. f|we present the temperature dependence of the drag scattering rate for a symmet- 
ric double quantum well system as in Fig. ^| but for different values of the matched electron 
density n (The second line from the top (n = 1.5 • 10 11 cm -2 ) is identical with solid line in 
Fig. |3|.). With increasing n the peak structure becomes broader while the maximum value 
decreases and the position T p shifts to higher temperature. This shift follows approximately 
the relation T p oc kp oc y/n (see the inset of Fig. |) which supports our findings that the peak 
position is related to the energy Tiskp which separates the regions of small (%u <C Tiskp) and 
large {Tiu ~ Tiskp) angle scattering. We have also studied the drag scattering rate as in Fig. 
|3] with fixed density n\ = 1.5 ■ 10 11 cm -2 but changing density n 2 at different temperatures. 
A peak occurs always at matched densities for which the momentum transfer is most effi- 
cient. It is seen for T = 6.5 K and most clearly pronounced at T = 2.1 K while for T = 1.1 
K it is only a small feature on an otherwise monotonous function of n 2 . This strikingly 
different behavior for T = 1.1 K as compared to T = 2.1 and 6.5 K which is similar to the 
experimental finding in Refs. |l3|Ji"4] , can be attributed to the changing dominance of small 



and large angle scattering with temperature and n 2 - 

In conclusion, we have calculated the drag scattering rate between spatially separated 
electron layers due to the virtual phonon exchange via DA and PA mechanisms We have 
shown that the frictional drag due to PA coupling is dominating. Assuming infinite lifetime of 
the intermediate virtual phonon states, the frictional drag does not depend on the interlayer 
separation. According to our calculations the temperature dependence of the drag scattering 
rate divided by T 2 exhibits a pronounced peak both at matched electron density of two 
layers and at temperatures where transition from small to large angle scattering occurs. 
These results are in agreement with the finding of J. Gramila, et al.0. 
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FIGURES 

FIG. 1. Schematic view of the double-layer system with the characterizing parameters (de- 
scribed in the text). 

FIG. 2. Diagram of the effective electron-electron interaction between layers 1 and 2 as medi- 
ated by exchange of virtual acoustic phonons. 

FIG. 3. Drag rate divided by the squared absolute temperature versus temperature calculated 
for the symmetric GaAs/AlGaAs double layer system of Ref. 13. The contributions of deformation 
potential (DA, dotted line) and piezoelectric coupling (PA, dashed line) are shown together with 
the sum (solid line). 

FIG. 4. Drag rate (sum of DA and PA) divided by the squared temperature as a function of 
temperature plotted for the symmetric GaAs/AlGaAs double-layer system as in Fig. 3 but with 
different electron concentrations. Inset shows the peak positions (squares) corresponding to drag 
rate data from this figure as a function of square root of the matched electron densities. Solid line 
is a guide to the eye. 

FIG. 5. Drag rate (sum of DA and PA) divided by the squared temperature for fixed first layer 
electron density n\ as a function of the second layer electron density ri2 at different temperatures. 
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